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Abstract. In this paper the regularity of optimal transportation potentials de- 
fined on round spheres is investigated. Specifically, this research generalises the 
calculations done by Loeper, where he showed that the strong (A3) condition of 
Trudinger and Wang is satisfied on the round sphere, when the cost-function is 
the geodesic distance squared. In order to generalise Loeper's calculation to a 
broader class of cost-functions, the (A3) condition is reformulated via a stere- 
ographic projection that maps charts of the sphere into Euclidean space. This 
reformulation subsequently allows one to verify the (A3) condition for any case 
where the cost-fuction of the associated optimal transportation problem can be 
expressed as a function of the geodesic distance between points on a round 
sphere. With this, several examples of such cost-functions are then analysed to 
see whether or not they satisfy this (A3) condition. 



1. Introduction 

In the optimal transportation problem, one seeks to find an optimal mapping from 
one mass distribution to another, such that a given cost functional is minimised over 
all mappings which satisfy a change of variables criterion. Formally, one is given 
two topological spaces Q~, Q + , a cost-function c : Q~ x Q+ — > R and two radon 
measures [ir and on Or and Q+ respectively, with the goal of finding an optimal 
map, T : — > that minimises the functional 

/ c(x,T(x))dQ- (1.1) 
Jq- 

over all mappings which satisfy the change of variables formula: 

f (hoT)dfi~=f hdfi + , V/7 e C°(Q + ;R). (1.2) 

jU + is referred to as the push-forward of [ir by the map T, which is notationally 
denoted as T#/i~ = The optimal transportation problem was first proposed by 
Monge in 1781 | Mon81 | and, in the above formulation, poses many obstacles to being 
analysed via standard methods in the direct calculus of variations (see | Urb98l lvN08| 
for details). 

For the rest of the paper, it will be assumed that dpr = f dx and d/j, + = gdy. 
Also, for the time being, it will be assumed that both Q~ and Q + are bounded sub- 
sets of K". Even with these assumptions, none of the issues with Monge's original 
formulation become resolved. Indeed, it was not until Kantorovich | Kan421 [Kan48 | 
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introduced the dual formulation of Monge's original question in 1942, that any signif- 
icant progress on this problem was made. Through Kantorovich's dual formulation it 
can be shown that an elliptic solution of the boundary value problem: 

Det [c xx (x, T u {x)) - D 2 u{x)} = |Det [c xy (x, T u (x))]j ,^ ( *\ v x 6 ft", (1.3a) 

T U (Q-) = Q + . (1.3b) 

yields a potential of Monge's original transportation problem, in a sense that if T u is 
subsequently defined as solving 

Du = c x (x, T (x)) (1.4) 

for a.e. x s Q~, then T u is indeed a minimiser of ( |1.1[ ) satisfying the change of 
variables criterion embodied in (1.2) (See EvaOTl I\7il03l lvN08"1 and references therein 



for details.) ( 1.3a ) is an elliptic partial differential equation of a Monge-Ampere type 



and is called the Optimal Transportation Equation. The boundary condition in (1.3b) 



is often referred to as a natural boundary condition 

Several conditions need to be placed on the cost-function in order to ensure exis- 
tence and uniqueness of optimal maps. Specifically, denoting U 2 Q~ x Q + , if the 
cost-function c is real-valued in C A {U) and satisfies 

(Al) For any (x,y) e U and (p, q) e D x c(U) x D y c(U), there exists a unique 

Y = Y(x,p), X = X(q,y), such that c x (x,V) = p, c y (X,y) = q. 
(A2) For any (x,y) 6 f, 

Det [D 2 xy c] ? 0, 

where D\ y c is the matrix whose elements at the / th row and j th column is 

dxjdyj ' 

then existence and uniqueness (up to a constant) of optimal maps can be readily 
inferred from the proofs in [Caf96, GM95]. The regularity of optimal maps gener- 
ated by solutions of the Optimal Transportation Equation for general cost-functions 
has only recently come into full understanding in the Euclidean space setting with 
the research presented in | MTW05 |, | TW06 | , | l_oeara | and | Loearb |. Essentially, the 
cost-function, in addition to having (Al) and (A2) hold, needs to satisfy: 

(A3) There exists a constant Q > such that for any (x,y) e U, and £, t/ £ R" 
with £ _L 77 such that 

(c«' r c ilq c r , st - c ljM )c 5M c tJ U jVkVl < -C |£| 2 |r?| 2 , (1.5) 
where c,j(x,y) = 9 g : }g' y v ' > , and [c' J ] is the inverse matrix of [c,j]. 



in order to ensure that solutions of (1.3a) are locally smooth (provided f and g are 



smooth) |MTW05| . The analogous global regularity result presented in [TW06] only 
requires a degenerate form of the (A3) condition to hold that corresponds to when 
C = 0. This degenerate version of the (A3) condition is labelled (A3w) in accordance 
with the notation in | TW06 |. While the research in |MTW05irTW06 | showed that the 
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(A3) condition was sufficient for higher regularity of optimal transportation potentials, 
it was Loeper in | l_oeara| that showed that the (A3) condition was indeed necessary for 
these regularity results to hold. Specifically, Loeper showed that if (A3w) is violated, 
one can build a pair of C°°, strictly positive measures, supported on sets with the usual 
smoothness and convexity assumptions, so that the optimal potential is not even C 1 ; 
and thus, the corresponding optimal map is discontinuous. Before this work, it was 
not known whether the (A3w) condition was truly fundamental for potential function 
regularity or if it was simply a technical condition to make the a priori estimates in 
1MTW05I and 1TW06] work. 

Remark. It should be noted that in addition to the above criterion on the cost-fuction, 
additional convexity and smoothness conditions need to be satisfied by dQ. ± in order for 
the global regularity estimates in [TW06] to carry through; this will not be elaborated 
upon in this paper. 

The author would like to thank Neil Trudinger and Xu-Jia Wang for their discussions 
and advisement regarding the research contained in this paper. 

2. Some Properties of the (A3) Condition 

In order to carry out the stereographic reformulation later in the paper, some ele- 
mentary properties of (A3) will be need. Specifically, the (A3) condition is symmetric 
in the x and y arguments, in addition to being invariant under coordinate transforma- 
tions. For clarity, these properties are derived below. 



To start off, it will be shown that the (A3) condition is symmetric in x and y. To 
demonstrate this, £ is defined as 



Using this definition to rewriting (1.5) with £ replaced by £, one sees that 

(c™Cij, p C q , rs - C U , rs )d- t cl- h C r - k d i - l Z t t h 7lkVl > Q|£| 2 M 2 , 

with a modified orthogonality criterion of 



0. 



(2.1) 



(2.2) 



The symmetry of x and y in the (A3) criterion is now evident in the rewritten form 
embodied in flzTl) and (pT2|). 



Remark. The symmetry of |2. 1| only makes sense when the mixed-Hessian of c exists 
and will still hold if the constant Q degenerates. This will be important in the forth- 
coming calculations, where costs that depend on geodesic distance between points on 
a round sphere will be considered. Specifically, these costs have an ill-defined mixed- 
Hessian for argument pairs that are antipodal points of each other. This situation and 
its consequences on the results of this paper will be elaborated upon in Section [5] 



Next, the expression depicted in (1.5) will be reduced to a less-cumbersome form, 
that also gives a better intuition as to what the (A3) condition itself actually means. 
Calculating, one see that 

Dp k c,j{x,Y(x, p)) = c u , q 



■ D Y q 

u Pk ' 



(2.3) 
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where the definition of Y (as stated in (Al)) has been utilised to gain the second 
equality. Differentiation of (2.3) subsequently yields 

D 2 piPkCu (x, Y(x, p)) = {c ihqr c qM + Cij - q cf) c rJ . (2.4) 

To proceed, a reduction of the term c q - k is required. From previous notational defi- 
nitions, it is understood that 

q, q c"' k = df. 

Differentiating this relation immediately produces the following relations: 

c' k J {x,y) = D Xk c'J{x,y) 

= c rj c kq Ax, y), 
c[i{x,y) = D yk c i Ax,y) 

= -d-"c rJ c q<kr (x.y). 



Combining ( |2.4[ ) with ( |2.5| ), it is now observed that 
D PlPk Cjj(x,Y(x,p)) 



-u.q 



Cr.st — Qj,st)c 



Thus, the (A3) condition is equivalent to 

D pl p k Cij(x,y)^jVkVi < -C |£| 2 M 2 . 
for a positive constant Q. 



(2.5) 



(2.6) 



With (2.6 ), it is now a straight-forward calculation to verify that the (A3) condition 
is also invariant under coordinate transformations. Fixing y, consider an arbitrary 
change of coordinates in x given by 

g{x) = x'. 
An elementary calculation shows that 

c l (x,y) = [D l g q ]c q (x',y) 

c u {x,y) = {[D ig ^[D jg r ]c qr {x',y) + [D ug s }c 5 (x' , y)) (2.7) 



From (Al) and (2.7), it is observed that 

p' k = c k {x',y) 

= [D k gT 1 c q (x,y). 



Using this in (2.7) yields 

c (J (x',y) = (lD i g*][Djg r ]c qr (x',y) + [D,^) 



(2.8) 



(2.9) 



From here, one can use the chain-rule along with the relation in (2.8) to deduce 

D Pk = [D k g q ]~ 1 Dp* (2.10) 



Thus, from (2.9) and (2.10), it has been shown that the left-hand side of (2.6) is 
transformed to 

D P(Pk c, 7 (x,y)[D,^][D J /][D,<g s ]- 1 [D,g t ]- 1 ^sr ?t . 

Redefining £, as [D,g q ]i q and rj, as [D l g q ]~ 1 r} q , finally shows that the (A3) condition 
is invariant under any change of coordinates, as the orthogonality criterion £ _L 77 is 
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preserved. 



3. Motivation 

In addition to showing the (A3w) condition to be sharp, in | LoearalH-oearb | Loeper 
also studies the correlation between curvature and the Optimal Transportation Prob- 
lem on Riemannian manifolds when the cost-function was of the form 

d 2 (x,y) 

c(x,y) = — (3.1) 

where d(x,y) represents the geodesic distance between points x,y e M". One of 
the interesting aspects of this research was that Loeper verified the theory therein 
(in the case where M" was taken to be a round sphere) by explicitly calculating that 



the cost-function depicted by (3.1) satisfies the non-degenerate form of the (A3) 



condition on round spheres. Specifically, it was shown that 

Dl Pk c l j{x,y)iiijT}m^-\, asx^y, (x, y) e Q~ x Q+ C § n . (3.2) 



This is particularly interesting as the cost-function in (3.1) only satisfies the (A3w) 
condition in the Euclidean case but not the stronger (A3) condition (see |MTW05l 
ITW06n . 



The disadvantage of Loeper's calculation in | Loearbl [Loeara | is that it relies heav- 
ily upon the specific geometric properties of the specific cost-function c(x,y) = 
\d 2 (x,y). In this paper, this calculation will be generalised on §" to include cost- 
functions which are arbitrary functions of the geodesic distance with respect to the 
constant curvature metric on §". 



4. Main Results 

The main results of this paper centre around the verification of the (A3) criterion 
for various cost-functions (on round spheres) having the general form 

c(x,y) = f(d(x,y)), (4.1) 

where d(x,y) represents the geodesic distance between x and y on the round sphere. 
The most important example thus considered is when c(x,y) = |c/ 2 (x,y). In this 
particular case, the following theorem (originally proven by Loeper and published in 
| Loearb | and subsequently improved upon in |KM08l [FRarl IFRVOQ I) is verified: 

Theorem 4.1 f |Loearbl IKM08I [FRarl IFRV09| ). Given S" be an embedded sphere 
in R n+1 with arbitrary radius, equipped with the round metric with an associated 
Riemannian geodesic distance d, if c{x,y) = |d 2 (x,y), then one has that c satisfies 
the strong (A3) condition on § n+1 x § n+1 \ {(x, -x) : x e §"}. 

The results in this paper extend the calculations in [Loearbj . Indeed, the precise 
constant for which the (A3) condition is satisfied, for the situation where c(x,y) = 
|c/ 2 (x,y), will be calculated for round spheres of arbitrary radius. Loeper calculated 
this constant for the case where R = 1; and the results in this paper indeed verify this. 
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Closely related to c(x,y) = |c/ 2 (x,y), is the cost-function of 
c(x,y) = 2*W(%^). 

This scenario corresponds to the situation where d is taken to be the chordal distance 
between two points in §" +1 relative to the ambient Euclidean space, into which the 
round sphere is embedded. For this case, the following analogous result to Theorem 
|4. 1| will be proven: 

Theorem 4.2. Given S" an embedded sphere in K" +1 with arbitrary radius, equipped 
with the round metric with an associated Riemannian geodesic distance d, ifc(x, y) = 
2R 2 sin 2 (%f^) , then one has that c satisfies the strong (A3) condition on 8" x §" \ 
{(x,-x)|xe §"}. 

In addition to the two above situations, a few other examples of cost-functions are 
studied later in this paper. 

It should be noted at this point that a cost-function satisfying the (A3) condition 
is not generally enough to guarantee higher regularity for Optimal Transportation 
potentials. Indeed, the Optimal Transportation Equation becomes singular for y e 
cut(x), for any cost-function depending on the geodesic distance between points x 
and y (specifically, | Det[D 2 y c]| blows up). Thus, in order to have potential function 
regularity, bounds on transport vectors, that prevent optimal transport maps from 
mapping points to their cut-locus, are also needed in addition to the (A3) verification. 
This geometric criterion will have implications in the forthcoming analysis, which will 
be discussed in Section [H 

5. Analysis of the (A3) Condition on §" 

The following formulation relies on several key observations and simplifications 
that are unique in the specific case of analysing the (A3) condition on round spheres. 
First, the derivation of the Optimal Transportation Problem (and hence the Optimal 
Transportation Equation) only depends on the measure-space structure associated 
with ft~ and Indeed, the only way geometry can come into the formulation of 
the Optimal Transportation Problem is if the cost-function is defined as having some 
explicit geometric dependence. The cost-functions analysed in this paper all depend 
on the geodesic distance between points on equipped with a round Riemannian 
metric. Thus, the underlying goal of the following calculation is to derive an explicit 
expression for the geodesic distance between two arbitrary, fixed points x and y lying 
on a round sphere. As will be discussed later, the round sphere is one of the very 
few manifolds (to the author's knowledge) where the geodesic distance between two 
arbitrary points can be explicitly represented. The rest of the simplifying observations 
will be stated as needed in the following formulation. 

Ultimately, the analysis of the (A3) condition will be carried out in R" (that is, 
in a local chart) with an associated, modified cost-function that yields an equivalent 
Optimal Transportation Problem, as compared to the one originally defined on As 
mentioned earlier, it can not be assumed that Optimal Transportation maps do not 
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move points to their cut-locus. Thus, in order to analyse the Optimal Transportation 
Problem on a manifold, via a single local chart, requires another geometric criterion or 
observation. If bounds on transport vectors exists such that optimal maps stay away 
from the cut-locus of particular points, then no additional criterion need be placed on 
the problem. As a broad class of cost- functions will be analysed in this paper, no such 
bounds have been proven for all the examples thus contained. Therefore, it will be 
assumed that cut(ft~)nQ + = 0, in addition to the standard Optimal Transportation 
hypotheses placed on Q~ and unless stated otherwise. 



Since only cost-functions will be considered that have the form depicted in (4.1), 
it is possible to analyse the Optimal Transportation Problem in a local chart with a 
modified cost-function whose associated Optimal Transportation Problem in Euclidean 



space is equivalent to that of (4.1) on S". Such a local chart and modified cost- 
function are formulated via a stereographic projection of §" embedded in R" +1 onto 
the some arbitrary tangent space of With this projection, an explicit expression for 
d(x,y) can be derived in terms of x and y: the projected coordinates on an arbitrary 
tangent space of S". 

5.1. Stereographic Projection. Given that it is assumed that cut(!T2~) n Q + = 0, 
there will be no geometric issues with the forthcoming calculations being carried out in 
a single chart of The possibility of relaxing this criterion will be discussed at the end 
of this paper, in Section [7] As the following stereographic formulation is tantamount 
to analysing the (A3) condition through a specific coordinate transformation, it is 
recalled that the (A3) condition is invariant under coordinate transformation according 
to the calculations in Subsection [2] and thus, the result of the following calculations 
will hold in general. 

5.1.1. The Half-Sphere Stereographic Projection. To begin, the modified form of 
( |4.1[ ) is first derived using the stereographic projection on the half-sphere depicted in 
Figure [T] Following this derivation, it will then be described how half-sphere projec- 
tion can be modified to project all points on the sphere to the plane, excluding the 
antipodal point of the intersection of the projective plane and sphere. 

Utilising the ambient Euclidean geometry of E" +1 , it is an elementary calculation 
that yields 




R 2 + x ■ y 



VR 2 + |*| V* 2 + \9\' 



(5.1) 



where the origin of R n+1 is set to coincide with the centre of the sphere without any 
loss of generality. 



From this calculation, the analysis of (4.1 ) on 8" is reduced to studying the Optimal 



Transportation Problem associated with the cost-function 

c(x,y) = f ( /?-arccos ( . /? +* "-> ? ]] ( 5 2 ) 

V W R2 + IWR 2 + lyl 2 ;; 
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Figure 1. Stereographic projection for the half-sphere 



between fi~ and fi+ on the local chart; and thus, (5.2) will now be analysed in the 
Euclidean setting. 



Remarks. 



(1) Even though the ambient Euclidean geometry of W +1 is used to derive (5.1), 
having the sphere embedded in R n+1 is not technically required nor is the 
association of the local chart with a particular tanget plane of §". Indeed, 
these notions are intuitive conveniences, in the sense that the equivalence 
of the Euclidean formulation of the Optimal Transportation Problem, to the 
original problem on is immediate by this explicit stereographic projection. 
A key point to the above calculation is that there is a simple relationship 
between 6 and the geodesic distance between two arbitrary points on §". If 
one were to project rays from a point on S" instead of its centre (as in the case 
of a full-sphere stereographic projection), there would be no simple relation 
between the angle of the projected rays and the associated geodesic distance 



(2) 



between two arbitrary points. The strength of (5.1) is that it is independent 



of the particular choice of tangent plane on which the stereographic projection 
is performed. 

5.1.2. The Full-Sphere Stereographic Projection. In order to analyse the (A3) condi- 
tion between any points on a round sphere (excluding pairs of antipodal points), the 
half-sphere stereographic projection must be slightly modified. Specifically, the coor- 
dinate transform itself has to be changed such that the projected point is associated 
with a ray emitted from the centre of the sphere but subsequently has been refracted 
back toward the tangency point at the sphere's surface. This projection is depicted 
in Figure [2] where x is taken to x* through this new coordinate transformation. 
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Figure 2. Stereographic projection for the full-sphere. Note that 
the projection ray corresponding to x has only been included for 
comparison to the half-sphere projection. 



Using basic trigonometry, it is readily calculated that 

|x*| = 2R - tan 

loreover, one has that 

x = | 1 - tan" 



Combining the above two relations, yields 



AR 2 



4R 2 - |x*| 2 



x . 



This final relation allows one to rewrite (|5.2|) as 
c ( x *.y*) = f \ R ■ arccos 



1 + (4R 2 -|x'j L2 ')(4R 2 -|y| 2 )( X * ' y*) 



(5.3) 



(5.4) 



(5.5) 



(5. 



The advantage of (5.6) is that the relation holds between any two points on a 



sphere, excluding antipodal point pairs, for any arbitrarily chosen projective tangent 



plane. It is important to note that the construction (5.6) from the half-sphere stere 



ographic projection is only done for expositional clarity; indeed, one could derive (5.6) 
with no use of the previous half-sphere projection. 

Notation. The forthcoming calculations will be based on the full-sphere relation in 



(5.6 ), but x. and y will be used instead x* and y* to indicate the full-sphere projected 



coordinates for the rest of the paper. 
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5.1.3. Stereographic Reformulation of the (A3) Condition. With the relation (5.6), 
it is now possible to formulate a new expression for the (A3) criterion in the current 
scenario. To do this, the geodesic distance on 8" is first differentiated with respect 
to the projected Euclidean coordinate. Subsequently, one may take the resulting 
expresion and set x to be the origin of the local chart without any loss of generality. 
This differentiation and origin selection ultimately give 



4,(x,y) 



-y> 
\9\ 



(5.7) 



Differentiating twice with respect to the projected, Euclidean coordinates and subse- 
quently choosing x to coincide with the origin yields 

\9Y 



4,x,(x,y) 



' 1013 



(5. 



Using these two relations, one can now readily differentiate (5.6) twice to find 



(5.9) 



li>l 2 V \9\ J ' IJ \9\ ' 

where x = has been assumed without loss of generality (since the choice was made 
after differentiation). 



To proceed, both y and d need to be represented in terms of the transportation 
vector, p, which is defined in the formulation of the Optimal Transportation Problem 
by 

p := Vxc(x.y), 



where c(x,y) is recalled from (5.6). Note that the vector notation on p will be 
suppressed from this point forward without loss of clarity. From this expression and 



(5.7), it is ascertained that 



Pi 



-f'(d) 



y, 



that is, 



\f'(d)\, and 



\9Y 

El 
\p\ 



y, 



(5.10) 



On the other hand, it is clear from Figure [T] and (5.3), with x = being assumed, 
that one has 

2R-tan (J^j = |y|. (5.11) 

Moreover, since cut(Q~) n ft + = is assumed, it is guaranteed that ^ < 7r. 

Remark. Again, it is important to note that the selection of the projective tangent 
plane is chosen and fixed only after the x derivatives on the cost function are per- 
formed. Once the tangent plane is fixed, x can then be chosen as the origin in that 
plane. 



Using both (5.10) and (5.11) , one may rewrite (|5.8[) as 



PiPj 



(f"(d)-fi(d)-£(d))+5 u fi(d)-£(d), 



(5.12) 
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where 



£(d;R) :-- 



2R-tan 



2RJ 



Defining the origin of the local chart to coincide with x effectively shifts all p depen- 
dence onto the y; and thus, differentiations with respect to p variables may subse- 
quently be applied directly to (|5.12h with the understanding x = 0. 



From (5.10) it is seen that d does indeed have p dependence. To differentiate d 



with respect to p, one must first rewrite the expression in (5.10) as 

d (sgn(f')|p|). 

Differentiating implicitly, it is seen that 

_ sgn(f')p, 

P '~ f"\p\ 
Pi 

~ Ji . f» ■ 



(5.13) 



where (5.10) has been used directly to ascertain the second equality, and the argu- 



ment on f has been suppressed without loss of clarity. 



Now, the final form of the (A3) condition (associated with the general cost-function 
f(d(x,y))) is able to be written. Taking two arbitrary, unit vectors £ and r\ such that 
£ _L 77, one has that 



Di 



(p-v) 2 


(£ 


IpI 2 




f" 




JTfii ~ 


friz 



£_ 

T 7 



£"f 



IpI 2 

£'f'f" 



f" 
flfT, 



f" 



+ 8 



f" 



fin 
£"f 



f" 3 



£' £ 

• 1 — 

fii fn 1 fi 

(p-ifip-nf (f»»_ 

f//2 



IpI 4 



f/2 fia 



£/ £lflfm 



(5.14) 



Remark. (5.14) gives the explicit, radial-scale dependency of the (A3) term in the 
current scenario. This has the interesting implication that it is possible to design 
cost- functions of the geodesic distance corresponding to the round sphere, so that 
for certain radii, the cost is strictly (A3) and for radii, it is not (A3) at all. This radial 
scale dependence is demonstrated explicitly in some of the forthcoming examples in 
Section [6] 



For convenience, (5.14) is used to define four functions of d: Pi(d), Piid), P$(d) 
and P${d), such that one may write 

(p-0 2 



Dlp^xjid^jVkW 



Pi(d) 



IpI 2 



-P 2 (d) 



ipi : 



IpI 4 



(5.15) 



To analyse (5.15), various orientations of the vector p relative to £ and 77 must now 
be considered. These calculations can effectively be reduced down to analysing four 
cases. 

Case I:: If p ± span (£, 77), then 



D 



PlPk 



c mj (dMjVkVi = Pi(d) =: Oi(d). 



(5.16a) 
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Case II:: If p|| £, then 

D 2 PiPk c XiXj (d)UjV k W = Pi(d) + P 2 (d)=:0 2 (d). 
Case III:: If p || 77, then 

D 2 plPk c m (d)^jV k W = Pi(d) + P 3 (d)=:0 3 (d). 
Case IV:: If p e span (£, 77) and p ■ $ = p • 77, then 

P 2 (d) , P 3 (d) , P 4 (d) 



(5.16b) 
(5.16c) 



£?2 



,cx,&j(d)^ijrj k rji = Pi(d) 



4 (d). (5.16d) 



2 2 4 

Here, four new functions of d are again defined for clarity. Subsequently, the combined 
negativity of 0\, O2, O3 and O4 is tantamount to the (A3) condition being satisfied. 
This is formally stated in the following lemma: 

Lemma 5.1. Given a cost-function of the form c(x,y) = f(d(x, y)) and S" an 
embedded sphere in R n+1 with arbitrary radius R, equipped with the round Riemannian 
metric with an associated geodesic distance d, then the following statements are true: 

(i) IfOj(d) < C < for i E {2,3,4} for all d E (0, Rir) and n > 2, then the 
strong (A3) condition is satisfied with constant C, for the cost-function. 

(ii) IfOi(d) < for i E {2, 3, 4} for all d E (0, Rir) and n > 2, then the (A3w) 
condition is satisfied for the cost-function. 

(iii) lfO,{d) <C < for i E {1, 2,3,4} for all d E (0, Rtt) and n>2, then the 
strong (A3) condition is satisfied with constant C for the cost-function. 

(iv) IfOj(d) < fori e {1,2,3,4} for all d e (0, Rtt) and n > 2, then the (A3w) 
condition is satisfied for the cost-function for n > 2. 



With the above calculations and recalling the statements made after (5.11), veri 



fication of the (A3w) condition is now reduced to verifying that (5.16a )-(5.16d ) are 
non-positive for d 6 (0, Rtt). Correspondingly, (5.16a )-(5.16d ) being strictly nega- 



tive indicate that the strong (A3) condition is satisfied. The primary benefit gained 



by the use of (5.16a)-(5.16d), is that one may uses these relations to easily generate 



explicit, analytic expressions that allow for straight-forward verification of the (A3) 
criterion for a general class of f. Such calculations will be the focus of the next 
section. 



6. Examples 



In this section, various examples of cost-functions will be analysed using (5.16a) 



(5.16d) to see if they do indeed satisfy the (A3w) or the (A3) condition. These 
examples include the some of the cost-functions analysed in | TW06 | and |MTW05| 
for the Euclidean case and are encompassed by the general cost-function 

c(x,y) = f(d(x,y)), (6.1) 

where d(x,y) represents the geodesic distance between x and y with respect to the 
underlying Riemannian manifold. 

Remark. The forthcoming results will be conveyed with the understanding that they 
hold for all d E (0, Rtt), unless otherwise indicated. 
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6.1. c(x,y) 



(x,y). The first example consider is the cost-function that has 
been the most studied out of all the possibilities encompassed by (6.1): c(x,y) = 



|d 2 (x,y). Indeed, this is the only cost-function for which the stereographic formu- 
lation may be applied without any extraneous geometric conditions, that serve to 
validate analysing the Optimal Transportation Problem in only one chart. This is due 
to the gradient estimate first proved in | McC01 | for compact Riemannian manifolds, 
which was then improved upon in the case of round sphere in [DL06]. To the au- 
thor's knowledge, there exists no such gradient bounds for the general cost-function 



depicted in (6.1). This estimate will now be briefly reviewed, as it will also motivate 



points of the discussion at the end of this paper. 

As was stated above, the scenario where c(x,y) = |d 2 (x,y) has been studied 
previously in | McC01 | in the context of Optimal Transportation on general Riemannian 
manifolds; and it has been specifically analysed on the round sphere in both [LoearbJ 
and | DI_06 |. In | McC01 |, McCann studies gradient mappings defined on Riemannian 
manifolds, which are mappings of the form 

G (m) := exp m (V m 0), 

where (p is the associated gradient-potential of the mapping. McCann shows that such 
maps are indeed optimal in the transportation of measures on Riemannian manifolds 
for the cost- function c(x,y) = |d 2 (x,y). In addition to this, McCann also shows that 
if is c-convex, then the length of it's gradient can not exceed the diameter of the 
manifold. Translating into the context of stereographic projections of round spheres, 



this means that the transport vector defined in (5.9) is not necessarily bounded in the 



projected local chart, but it is well-defined there. Delanoe and Loeper improved this 
result on S" by proving the following gradient bound in | DI_06| : 

Theorem 6.1 ( |DL06j ). Given <p '■ §" >-> K 3 c-convex function, such that 



(/7oG ) dVol = / (h-p) dVol 



(6.2) 

§", dVol) with any h e C°(§"), where dVol stands for the canonical 



for some p e L 

Lebesgue measure on S", then the following estimate holds a.e. 



\d<t>\ < 7T 




nVol(S") 



I 2 



2Vol(S"- 1 ) 



l/n 



(6.3) 



Remark. Recalling that the notation remarked upon in Section [T] the expression in 
(6.2) is often denoted as 

dVol = pdVol . 



Theorem 6.1 correlates to the transportation vector p, as defined in (5.9), being 



strictly bounded in the stereographically projected local chart. Thus, the assumption 
that cut(ft~) n n + = is not needed in the case where c(x,y) = \d 2 (x,y). The 
possibility of extending (6.3) to the general case of c(x,y) = f(d(x,y)), will be dis- 



cussed at the end of this paper. The (A3) condition will now be analysed for the case 
when c(x,y) = |c/ 2 (x,y), using the results of Section [5j 
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Using ( 5.14[ ), it is calculated that 
cos ((2/?)-^) 



Pi(d) 
P 2 (d) 
P 3 (d) 

P A (d) 



2Rc/sin((2R)- 1 d) 4R 2 sin 2 ((2/?)-^) 

_2 cos((2R)- 1 c/) 1 

"d 2 + 2/?dsin((2/?)- 1 d) + 4P 2 sin 2 ((2P)- 1 d) 



cos 



1 dcos((2R)- 1 c/) 
2Rc/sin((2R)- 1 c/) 4R 2 sin 2 ((2R)~ l d) + 2 8/? 3 sin 3 ((2R)^d) 



((2/?)- 1 d) 



(6.4a) 
(6.4b) 

(6.4c) 



8 dcos((2R)- 1 c/) 



cos ((2/?)- 1 d) 



e/ 2 8R 3 sin 3 ((2R)" 1 d) 2Rdsin ((2R)- 1 c/) 4R 2 sin 2 ((2/?)- 1 d) ' 

(6.4d) 



As stated in Subsection 5.1 (6.4a )-(6.4d ) are considered for d £ (0, Rtt). First, the 
limits as d — > are calculated to be 
2 



lim Pi (of) 

d-s-0 



limP 2 (d) = 0, lim P 3 (d) = 0, lim PJd) = 0. 

3P 2 d^O d^O d^O 



Using these limits to analyse Oi(d), 02(d), 03(d) and 04(d), it follows that 

2 



^ D l Pk c ^(d)UjVkVi - 3/?2 , 



This confirms the result of Loeper presented at the end of | Loearb| for the case where 
R = 1. Using ( |6.4a| )— ( [6~4d] ) , it is also calculated that 

o;(d)<0, Vc/G(0,R7r), VR>0, 

where / e {1,2,3,4}. Figure [3] below contains plots of both sets of functions P,(d) 
and 0;(d), in the case when R = 1. 




With the above calculations, it has been shown that the (A3) condition is indeed 
satisfied on the round sphere of any radius, for the cost of c(x,y) = ^d 2 (x,y). 
As noted earlier, this is in contrast to the Euclidean case where the cost-function 
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y| 2 only satisfies the (A3w) condition, as the (A3) term goes to as 



c(x,y) = \ 

y — > x (see | MTW05l ITWOG I). This contrast between the Euclidean and spherical 
cases indicates that the underlying geometry (manifested through the cost-function) 
does indeed affect the local regularity of the Optimal Transportation Problem. 



6.2. c(x, y) = 2R 2 sin 2 ( d ^R^ ) ■ If one were to consider the round sphere embedded 

in equipped with a Euclidean metric, the cost-function of 2R 2 sin 2 (j^d(x,y)) 
is equivalent to the example studied in Subsection|6.1|with d(x, y) taken to being the 



geodesic distance of the Euclidean space into which the sphere is embedded. This 
situation is depicted in Figure [4] on the next page. 




2Rsin (f) 



Figure 4. Chordal distance between points 



Again, using the results of Subsubsection [5. 1.3| it is calculated that 
1 



Pi(d) 
P 2 (d) 



4R 2 [l-2cos 2 ((2R)- 1 d)] 

4 sin 2 ((2R)- 1 c/) cos 2 ((2/?)- 1 d) 



(6.5a) 



AR 2 [1-8 cos 6 ((2/?)" 1 c/) + 12 cos 4 ((2/?)~ 1 d) - 6 cos 2 ((2R)- 1 d)] 

(6.5b) 



P 3 (d) = 
P A (d) 



(6.5c) 



2sin 4 ((2R)- 1 c/) cos 2 ((2/?)- 1 d) 



AR 2 [8cos 6 ((2R)- 1 d) - 12cos 4 ((2R)- 1 d) + 6 cos 2 ((2RY 1 d) - 1] ' 

(6.5d) 



From these relations, the following limits are calculated 
1 



lim PAd) 

d-s-0 



\imP 2 (d) = 0, limP 3 (d) = 0, lim PJd) = 0; 

d^O d^O d^O 
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thus, 



R 2 ' 



(6.6) 



Finally, an elementary calculation using ( 6.5a[ )-(6.5d ) in the definitions of 0-,(d) yields 
O;(d)<0 Vc/G(0,R7r) VR>0. (6.7) 



Given the properties depicted in (6.6) and (6.7), it follows that this cost-function 



does indeed satisfy the strong (A3) condition on the half-sphere. This is readily 
confirmed by the plots of P;{d) and 0,(d) in the case where R = 1, shown in Figure 
[U below. 



■ p„(d) 



\p 2 (d) 



\O t (d) 



Figure 5. P;{d) and 0;(d) corresponding to c(x,y) 
2/? 2 sin 2 with R = l 



6.3. c(x,y) = \fl — d 2 (x,y). This example gives a demonstration as to just how 
important R scaling can be in affecting the (A3) condition. This radial scaling de- 
pendence manifests itself through the limiting behaviour of the orientation terms as 
d — > 0. Using the general expression for 0,(d), one readily calculates that 

HmO,-((0=3^-l for / e {1,2,3,4}. 

Further calculation shows that for R > yfljZ, and diam(fT U < 1, the (A3) 
condition is satisfied. For R = \/2/3, the (A3) condition is satisfied if dist(fi~, Q. + ) > 
and diam(ft- U Q+) < 1. If one only has that dist(Q~, Q+) > and diam(r2" U 
Q + ) < 1, then only the (A3w) condition is satisfied. In other cases where R > a/2/3, 
the (A3) condition will be violated. 

If R < v/2/3 and diam(ft- U ft+) < 1 with dist(Q-,Q+) > h* , where h* > 
solves P%{h*) = 0, then c(x,y) = yl — d 2 satisfies the strong (A3) condition. This 
cost-function will violate the (A3) condition in all other cases where R < \/2/3. 
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6.4. c(x,y) = - d 2 (x,y). This example is closely related to the previous ex- 
ample, in that it also has a dependence on the radial scaling of the round sphere. In 
particular, one has that 

limO,(c/) = l-^2 for i e {1,2,3,4}. 

From this, it is observed that for R > \/2/3, the (A3) condition will be violated at 
least for some small values of d. Indeed, further calculation indicates that 

O;(d)>0 for ie {1,2,3}, Vd e [0,1], 

for any R > \/2/3, thus proving the (A3) condition will be violated in this case. 

If R < y/2/3 and diam(ft-UQ+) < m\n{h* , R, 1}, where h* > solves P 2 (h*) = 0, 
then the strong (A3) condition will be satisfied. In all other cases where R < a/2/3, 
the (A3) condition will be violated. 

6.5. c(x,y) = y/l + d 2 (x,y). This example bears some resemblance to the previous 
one; but in this scenario, there is no sign-symmetry broken by a variance in the radial 
scaling. Indeed, it is calculated that 

limO ( -(c/) = -r§2-l. for/ G{1,2,3,4}, 

along with 

O'(d)<0, i e {1,2,3,4}. 
In this scenario, the (A3) condition will be strongly satisfied for all values of R > 0. 

6.6. c(x,y) = — \/l + d 2 (x,y). In this example, it is calculated that 

O;(d)>0, ie {1,2, 3, 4}, 
for any R > 0. Thus, the (A3) condition will always be violated for this cost-function. 

6.7. c(x,y) = ±^d m (x,y). In this subsection, a generalisation of the example in 
Subsection 16.11 is considered: 



d m 

c( Xl y)=i-(x 1 y)m/0 1 c(x, y) = ± log(d(x, y)) m = 0. (6.8) 



Calculating as before, one finds the following 

(m - l)2Rsin ((2R)" 1 c/) cos ((2R)- 1 d) - d 
(m - l)4R 2 d m - 1 sin 2 ((2/?)" 1 c/) 



I/// 1 \ Jill I V I U LUJ U-l M LJ / u 

Oi(d) = ±- ^ — ; u — ^ (6.9a) 



m2Rsin((2R)- 1 d) -2dcos((2R)- 1 c/) 
° 2 ^ = T 2/^sin((2ft)-id) (6 - 9b) 



,m2/?sin((2/?)- 1 d) - 2d cos ((2R)- 1 d) 
(m - l) 2 8/? 3 d m " 2 sin 3 ((2/?)-^) 



Q 3 (d)=T - "J^.i_,_ : , 3 ' (6-9c) 



_ cos ((2/?)-^) cos ((2/?)-^) 

U 4(.oj - ±tt- ^7T7-^T^— 377^7-^T7; ± 



2(m - l) 2 8R 3 d m - 3 sin 3 ((2/?)- 1 d) 8Rc/ m_1 sin ((2/?)"^) 

6m -5 m 2 -2m + 2 (m - 2) 2 

4(m - l) 2 4R 2 d m - 2 sin 2 ((2/?)- 1 d) 2(m - l)d m " 4d 2 

(6.9d) 
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From (6.9a )-(6.9d ), the following series of observations can be made for various 
values of m and sign on the cost-function depicted in (6.8). 

6.7.1. (+), m < 0. In this scenario, the R scaling has no effect on the (A3) results. 
From (6.9a )-(6.9d ), the following is calculated: 

P/(d)>0, for ie {1,2}, 

Pj(d)<0, for j e {3,4}, 

\Pi(d)\ > \P 3 (d)l (6.10) 

for d e (0, Rtt). Analysing limits as y — > x, it is observed that 

0, for k e {1,2,3}, 



lim P k (d) 
lim PJd) 



-oo. 



(6.11) 



From (6.10) and (6.11), it is thus concluded 

Oi(d)>0, for /' e {1,2,3}, 

for d e (0, Rir); that is, the (A3) condition does not hold in this scenario. 

6.7.2. ( — ), m < 0. Although this case bears many similarities to the scenario de- 
picted in Subsubsection |6.7~l) it also has some notable differences. In particular, the 
radial scaling of the sphere does affect whether or not the (A3) condition is satisfied. 
However, results that are true for all R > will first be conveyed. Analysing the Pj(d) 
expressions, it is calculated that 

P/(o/)<0, for ie {1,2}, 

P 3 {d) > 0, 

\Pi(d)\ > \P 3 (d)\, (6.12) 
for d e (0, Rtt) and VR > 0. To proceed, the limits asy->x are calculated: 

0, for ke {1,2,3}, 

-oo. (6.13) 



lim P k {d) 

d->0 



lim PJd) 

d-s-0 



From (6.12) and (6.13), one has the following: 

O,(c/)<0, for ie {1,2,3}, 

for d e (0, Rtt). So far, no mention has been made of the 04(d) expression in the 
current analysis. This is due to the fact that this term can switch signs on the interval 
d e (0, Rtt), which is a result of Pn(d) becoming positive for some d > 0. (6.13) 
already shows that 04(d) — > —00 as y — > x; and thus, it is known that there exists 
some constant M such that d < M implies 4 (d) < by the continuity of O 4 (o'). 
Thus, the analysis of the 4 (gQ expression is reduced to analysing the equation 







-2gP cos 



d 

2R 



d 

2R 



(6m-5)2Rd 4 sin 
2(m - l)(m 2 - 2m + 2)8R 3 d 2 sin 3 
(m - l) 2 (m - 2) 2 8R 3 d m sin 3 (~ 



d 

2R 



l) 2 4R 2 d 3 sin : 



d 

2R 



cos 



(6.14) 
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the right-hand side is simply the numerator of (6.9d). For m < and d e (0, Rtt), 
the fourth term on the right-hand side of (16.14) is the only positive term. If we 



consider d e (0, 1] and m < 0, it is observed that 



2(m - l)(m 2 - 2m + 2)R 3 d 2 sin 3 



2R 



< 



(m- l) 2 (m- 2) 2 R 3 cTsin 3 



d 

2R 



(6. 



as 



,m-l) 2 (m-2) 2 > 2(m-l)(m 2 -2m + 2) for all values of m. Thus, by ( |6.13[ ) and 
Q6.14[) it is seen that 4 (d) < for c/ e [0, 1]. The actual value d* where 4 (d*) = 



can not be explicitly represented as ( 6.14| ) is a transcendental equation. This value 
of d* depends on both R and m. Even though an analytic representation of d* does 



not exist, it can be calculated from (6.14) that 



lim d* 

m— > — oo 



VR > 0. 



(6.16) 



This can be seen from (6.14) by noticing the fourth term on the right-hand side 
dominates all terms for large negative values of m and d > 1 and some R > 



fixed. However, (6.15) holds for all m < 0, R > and d < 1; thus, one can intuitively 



reconcile (6.16) from (6.15) without resorting to a limit calculation. (6.16) represents 



the manifestation of the R scaling dependency for this specific scenario. In particular, 
if R < i the (A3w) condition will be satisfied. Of course, this is not a sharp estimate; 
but such restrictions on the radius of the sphere will ensure the (A3w) condition is 
satisfied for any m < —2. Instead of restricting the radius of the sphere, an analogous 
restriction may be employed on the source and target domains: diam(fi~ U fi + ) < 1. 
This will also ensure that the (A3w) criterion is satisfied. In order for the strong (A3) 
condition to be satisfied, one of the aforementioned restrictions is required plus the 
criterion that dist(Q~,ft+) > 0, as 



limO,(ef) = 0, for ie {1,2,3}, 



which is a straight-forward consequence of (6.13) 



As was already mentioned, the limit (6.16) only corresponds to a sufficient condition 



on R for the (A3w) condition to be satisfied. The corresponding necessary condition 
on R can indeed be calculated numerically for given values of m; but such a condition 



can not be explicitly stated due to the transcendental nature of (6.14). However for 



a fixed R, this scenario undergoes a bifurcation in the m parameter whereupon the 
(A3) condition will be at least weakly satisfied for all R > 0. The value for m, that 
corresponds to this bifurcation, depends on R and is again only able to be calculated 
numerically. One may analyse this bifurcation point as R — > oo. Assuming R > and 
d > 1 is such that ^ << 1, the right-hand side of ( 6.14[ ) is approximated by 



(-2m 3 + 5m 2 - 4) d 5 



(d 7 ) 



By this approximation, it is seen that if m* < m < 0, where m* is defined as the 
the negative root of the polynomial —2m 3 + 5m 2 — 4, then for an arbitrary fixed 
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d e (0, R-k) 



> lim 



-(m- l^P/o^sin 



d 

2R 



cos 



R 



-2(m - l)(m 2 - 2m + 2)8P 3 d 2 sin 3 



2R 



It can be calculated that 



m 



-.7807764064, 



(6.17) 



which is a numerical limit for the (A3), m-parameter bifurcation as R tends to infinity. 
This is shown using elementary methods of calculus to conclude that for m > m* , the 



first three terms of (6.14) dominate the single positive fourth term on the right-hand 



side of (6.14) 



Remark. It indeed does require the sum of all three of these terms to dominate this one 
positive term in the 04(d) expression; thus, there is no simple asymptotic statement 



to justify this behaviour outside of analysing the sign of the right-hand side of (6.14) 
and it's derivatives. 



Again, (6.17) correlates to a limit of R tending to infinity. As R decreases, this 



bifurcation will happen for a value of m < m*. If R falls below the previously discussed 
bifurcation point dependent on m (we recall this has to be greater than then the 
(A3) condition is at least weakly satisfied. This combined with the previous statements 
regarding the bifurcation point in the R parameter fills out the current spectrum of 
results regarding this scenario. The results for this particular case are summarised as 
follows: 

• If R > ^ and m > m*, the (A3w) condition is satisfied for any Q~ and Q + 
such that cut(fT) n Q + = 0. 

• If R > i and m < m* , one must numerically check to see if (16. 141) has 



real zeros for d e (0, Rir). If this is the case, the smallest positive root of 



(6.14) represents an upper bound for the diameter of U Q + for the (A3w) 
condition to hold. 

• If R < \ the (A3w) condition will hold for any Q~ and fi + such that 

cut(fi-)nft+ =0. 

All these conclusions can be strengthened to having the (A3) condition satisfied, 
provided dist(Q",Q+) > 0. 



6.7.3. (+), m = 0. This, scenario is exactly the same as the case studied in Subsub- 
section |6.7.1| except that 



lim Pi (of) 



VR > 0, 



which is calculated from (6.9a). Thus, the (A3) condition will not be satisfied in this 
case. 
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6.7.4. ( — ), m = 0. Straight-forward calculations from (6.9a )-( 6. 9d ) indicate that 

O,-(cQ<0, for/ G{1,4}, 
Oj(d)<0, forje{2,3}, 

for d E (0, Rir). Indeed, analysing limits, it can be ascertained that that 

lim Oj{d) = -2, for iE {1,2,3}, 

d^O 

lim OJd) = —oo 

d^O 

and 

lim 0,(c/) = 0, for / E {2,3} 

d^R-K 

lim 0,(c/) < 0, for / E {1,4}. 

These results are true for all values of R > 0. Thus the (A3w) condition is satisfied 
independent of the radial scaling of the round sphere. To have the (A3) condition be 
satisfied, it is required that dist(cutQ~, > 0. 

6.7.5. (+), < m < 1. As with the cases where m < 0, the scenarios for when 
m > also have a complex structure in that both the varying of R and m have bifur- 
cations in all of the terms 0,(d) for / E {1,2,3,4}. The following exposition will be 
less explicit as compared to the m < cases, as the analysis is extremely similar to 
examples already presented in this section. 

In this case the (A3) condition will be violated for all values of R > 0, except in 
a very special situation. First, one can easily verify that 0\(d) > for d E [0, Rir); 
but there are values of < m < 1, R > and d such that Oj{d) < 0. Specifically, 
for m > m* and any value of R > 0, one has the (A3) condition being satisfied on 
§ 2 provided dist(Q~,Q+) > h* , where h* represents the second positive root of the 
equation 04(d) = 0. m* is the root of the equation 

lim 4 (Rtt) = 0, (6.18) 

and is approximately .806. If R ^ .071, then h* is determined by the smallest positive 
root of 02(d) = which is equivalent to the equation O^(d) = 0. 

Remarks. 

(1) It should be noted that 04(d) has a parameter dependence on m for cost- 
functions of the form c(x,y) = ±-^.d m (x, y). With this, it is clear what is 



meant by m* being the root of (6.18). 



(2) This is one of the special cases where the dimensionality of S" has an affect 
on the (A3) condition. Indeed, for n > 2, one must have Oi(d) less than 
zero, which is not the case in the current scenario. However, for n = 2 the 
Oi(d) is not considered, as it is clearly not possible to have a transport vector 
orthogonal to both arbitrary vectors 77 and £ with 97 _L £ in a two dimensional 
space. 
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6.7.6. (+), 1 < m < 2. In this case, 

O,(d)<0, for ie {1,4}; 
and 02(d) and 03(d) are monotone decreasing with 

lim 0,(d) = 00, for j e {2,3}. 

It can be calculated that for this scenario there exists h* e (0, Rtt) such that 02(h*) = 
03(h*) = 0. Thus, h* represents the minimal distance separation between fi~ and fi + 
for the strong (A3) condition to be satisfied in this particular case. This conclusion 
is independent of radial scaling. 

6.7.7. (+), 2 < m < 00. In this case, one has that 

Oj(d)<0, for ie {2,3}, 
and 02(d) and 03(d) are monotone decreasing with 

lim 0,(d) = 00, for j e {1,4}. 

The strong (A3) condition will be satisfied in this case if dist(Q~,fi + ) > h* , where 
h* is defined by the equation 4 (/?*) = (there is only one root of this equation in 
the interval (0, Rir)). It is a straight-forward calculation to see that h* e (0, Rir). 
The results of this scenario are invariant under radial scaling. 

6.7.8. ( — ), < m < 1. Unlike the case for m < 0, a change in sign is not necessarily 
tantamount to violation of the (A3) condition. Here one has that 

Oj(d)<0, for/ €{1,4} 

and 

Oj(0) < 0, O'j(d) > 0, for j e {2, 3}. 

It is calculated that there exists a h* e (0, Rtt) such that 2 (h*) = 3 (h*) = 0. 
Thus, if diam(fi- U Q+) < h* , then the (A3) condition holds. 

6.7.9. (-), 1 < m < 2. In this case, it is calculated that 

Oi{d) > 0, 

Oj(0) < 0, O'j(d) > 0, for j e {2, 3, 4}, 

for all d e (0, Rtt) and R > 0. Thus, the (A3) condition will not hold for this 
cost-function on §" for n > 2. On § 2 however, there exists h* e (0, Rtt) such that 
4 (h*) = 0. Moreover, if diam(Q~ U ft+) < h* , then the (A3) condition will be 
strongly satisfied on § 2 . 

6.7.10. (-), 2 < m < 00. Here the (A3) condition can not be satisfied for any R > 
with 2 < m < 00, as it is calculated that 

Oj(d)>0, for /e {2,3}, 

for all d e (0, Rtt). 
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Remarks. This concludes a particularly long example that demonstrates the use of 
Lemma [5~Tl The complex structure that correlates to the above set of scenarios is due 
to the fact that one must analyse four different expressions for various orientations 
of £ _L r\ relative to the transport vector p, while varying two different parameters: 
the exponent m and the radius of the round sphere R. It has been demonstrated 
in this example that bifurcations in the (A3) behaviour happen relative to both these 
parameters for different orientation vectors, thus resulting in a large group of scenarios 
and corresponding results. 



7. Conclusions 



The full generality of (4.1) has been able to be studied at the expense of utilis- 
ing an inherently non-intrinsic approach in the current set of calculations. However, 
the use of the stereographic projection as the centrepiece of the analysis presented 
in this paper is uniquely powerful in the context of round spheres for two main rea- 
sons. First, the stereographic formulation is rotationally invariant on indeed, the 
analysis in this paper was carried out on a arbitrary fixed x£§" against a variable 
Optimal Transportation target y. The second (and most important) simplification 
the stereographic formulation affords is the ability to explicitly represent the geodesic 
distance between two points on a sphere in terms of the projected coordinates; that 



is, (5.1 ) is valid as a representation of the geodesic distance between points x, y e §" 
in terms of the projected coordinates x, y E W. Such an explicit and technically man- 
ageable representation of geodesic distance on a general Riemannian manifold is rare 
in any coordinate system on may choose for a chart on that manifold. Indeed, to the 
author's knowledge, representations of geodesic distances on even an ellipsoid result 
in the necessity to use highly esoteric special functions based on implicit or integral 
representations. 

As the (A3) analysis in this paper reduces to the analysis of covariant derivatives of 
f(d(x,y)), nothing can be calculated without an explicit representation of geodesic 
distance. Thus, given the above comments, it is computationally difficult to extend 
the methods in this paper to other Riemannian manifolds beyond that of the round 
sphere. 



As stated in Subsection 6.1 verification of the (A3) criterion is only one part 
to proving the regularity of potential functions associated with certain costs. The 
other part of proving regularity lies in the existence of gradient estimates analogous 



to those presented in Theorem 6.1 for cost-functions other than ^d 2 (x,y). To 
circumvent this, the assumption that cut(Q~)nfi + = has been made throughout the 
paper. Without this assumption or a gradient estimate, the stereographic formulation 
becomes invalid as a point may be mapped to it's cut-locus on the sphere; and thus, 
move outside of the chart where the analysis was performed. In the case where f'(d) < 
0, one can not escape the assumption that cut(Q~) n Q + = in the stereographic 
formulation. In this situation, one only need to consider a case where cut(Q~) = fi + 
to observe that the optimal mapping correlates to a mapping that takes every point 
in x e and maps it to that point's particular cut-locus. Thus, to analyse cases 
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where f'{d) < free from the assumption that cut(ft ) n Q + = 0, requires the use 
of geometrically intrinsic methods. 
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